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Abstract 

We present a systematic approach to deriving normal forms and related amplitude equa- 
tions for flows and discrete dynamics on the center manifold of a dynamical system at local 
bifurcations and unfoldings of these. We derive a general, explicit recurrence relation that com- 
pletely determines the amplitude equation and the associated transformation from amplitudes 
to physical space. At any order, the relation provides explicit expressions for all the nonva- 
nishing coefficients of the amplitude equation together with straightforward linear equations 
for the coefficients of the transformation. The recurrence relation therefore provides all the 
machinery needed to solve a given physical problem in physical terms through an amplitude 
equation. The new result applies to any local bifurcation of a flow or map for which all the 
critical eigenvalues are semisimple (i.e. have Riesz index unity). The method is an efficient and 
rigorous alternative to more intuitive approaches in terms of multiple time scales. We illustrate 
the use of the method by deriving amplitude equations and associated transformations for the 
most common simple bifurcations in flows and iterated maps. The results are expressed in 
tables in a form that can be immediately applied to specific problems. 
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The behavior of physical, chemical, and biological systems near the onset of dy- 
namic instabilities has considerable importance. Here the dynamics can be described 
efficiently in a low-dimensional space of "amplitudes" by using so-called amplitude 
equations or normal forms. Amplitude equations like the complex Ginzburg-Landau 
equation have been widely used to describe the qualitative behavior of such systems. 
The paper provides a general method for setting up all the necessary equations to 
solve a given physical problem quantitatively in terms of amplitude equations. 

For example, chemical oscillations and waves can be described by differential equa- 
tions which near the onset of oscillation can be approximated in terms of complex 
amplitudes. This reduces the state space from the (possibly high-dimensional) con- 
centration space to a complex plane. Such reductions usually result in immense savings 
in computational effort. To describe the dynamics in terms of amplitudes one must 
first of all determine the coefficients of the amplitude equation. For a quantitative 
study and comparison with experiments, one must also determine the transformation 
from amplitudes to the actual concentrations of all the chemical species. The paper 
provides tables of formulas for calculating all relevant coefficients. For instabilities not 
covered by the tables, the formulas can be calculated by a simple recurrence relation, 
the main result of the paper. 

1 Introduction 

The present work grew out of a need to extend the complex Ginzburg-Landau equation for a chemical 
reaction-diffusion problem In the process it became clear that a systematic and efficient 

method of deriving amplitude equations and associated transformations for local bifurcations is 
lacking. We provide such a method for an important class of bifurcations in the present paper. 

Briefly, the idea of an amplitude equation for an oscillatory system is to represent the state 
of the system through a complex amplitude (real amplitude and phase) which varies much more 
slowly in space and time than the original dynamical variables. A slow variation of an amplitude 
as compared with much faster oscillations apply near a supercritical Hopf bifurcation. In fact, the 
concept of an amplitude can be generalized to non-oscillatory modes near local bifurcations. 

A large difference of characteristic times of variation of amplitudes versus period of oscillation, 
is the basis of traditional "methods of multiple times" J^-|^] to derive amplitude equations. Unfortu- 
nately, such methods are somewhat intuitive (as opposed to rigorous) and are not easy to generalize. 
The property utilized in amplitude equations is the critical slowing down of modes associated with 
a bifurcation. This means that amplitude equations can be related to normal forms describing the 
essential features of the motion of a state point in the center manifold at a bifurcation. In this paper, 
we suggest a precise and general definition of amplitude equations associated with local bifurcations 
of a flow or iterated map and develop a systematic method of deriving them. We emphasize that 
the present work does not require any use of normal form theory. In fact, the derivation of the 
traditional normal form theorem appears naturally as part of the theory developed in this paper. 

We restrict the discussion to bifurcations associated with spatially homogeneous systems and 
exclude e.g. Turing bifurcations. Consequently, we shall treat systems defined by ordinary differen- 
tial equations or iterated maps. The theory can be generalized to spatially inhomogeneous systems 
e.g. with a linear diffusion operator. 

The idea of the method proposed is the following. Consider a dynamical system at a local 
bifurcation. The motion of the state point in the r-dimensional center manifold W c is described 
by r (generalized) amplitudes Zj (j = 1, . . . , r). These are coordinates of a point z in the center 
subspace E c with respect to a basis of eigenvectors of the linearized vector field or map J. We 
assume that J has r linearly independent eigenvectors in E c , but make no assumptions about the 
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other modes. 

For a flow, we derive a set of differential equations for the amplitudes Zj, the amplitude equation, 
and a transformation h(z) from E c to W Q so that the image h(z(t)) describes the dynamical 
evolution of the state point on W c when z(t) is the appropriate solution to the amplitude equation. 

The transformation h(z) is determined by the condition that as many coefficients of the ampli- 
tude equation as possible must vanish, i.e. it is a special type of normal form transformation. The 
non- vanishing coefficients of the amplitude equation are determined together with the transforma- 
tion h(z), order by order in an iterative process. The method is very efficient — in one process, it 
determines the amplitude equation (in normal form) and its unfolding as well as the center man- 
ifold. Initially, all that is required is to determine the non-hyperbolic eigenvalues of the Jacobian 
J and the corresponding right and left eigenvectors. (Thus, it is not necessary to transform J into 
Jordan block form.) At each order of iteration, the only non-trivial calculations are the solution of 
a few systems of linear equations. The center manifold is then locally determined by the transfor- 
mation h(z) that immediately provides the physical solution from that of the amplitude equation. 
Other representations of the geometry of the center manifold can easily be obtained from h(z). The 
method for iterated maps (including Poincare maps associated with periodic solutions to flows) is 
similar to that just described for flows. 

The method reported here makes it straightforward to solve a given physical (chemical, biolog- 
ical, . . . ) problem in terms of an amplitude description. Once the general explicit formulas have 
been derived for the type of bifurcation in question, the coefficients of the amplitude equation and 
the transformation h(z) can be calculated for the particular problem directly in terms of the vector 
field (or map). So the amplitude equation for the particular physical problem can be written down 
immediately, and any solution to it in terms of amplitudes can then be expressed in terms of real 
physical quantities through the transformation h(z). To illustrate the method, we derive general 
formulas for the most common simple bifurcations for flows as well as iterated maps. A more chal- 
lenging application of the method to derive an amplitude equation for an extension of the complex 
Ginzburg-Landau equation in the context of fold-Hopf bifurcations in reaction-diffusion problems 
will be published separately. 

A method similar to the one developed here has been obtained previously by Bruno (and 
see the references therein) for the general case. However, the restriction to semisimple critical 
eigenvalues allows us to derive results for any such bifurcation in particularly simple form. These 
results, a recurrence relation ( p8|) with ( p9| ) for flows and (H) for iterated maps, provide explicit 
expressions for the amplitude equation (or normal form), fl31[) and ( |55| ) respectively, and explicit 
linear equations for the transformation from amplitudes to physical space, (30) and ( |56|) respectively. 
It is entirely straightforward to obtain all the nonvanishing coefficients of the amplitude equation 
and the associated transformation, order by order to any desired order, and to apply these to specific 
physical problems. 

Explicit results for a number of specific bifurcations, applicable at the bifurcation points, can 



also be found in the recent book by Kutznetsov 



. Other approaches are described for example by 



Guckenheimer & Holmes |9(|, Coullet & Spiegel jl0|, Elphick et al. and by Golubitsky et al. |12| . 
For a more physically oriented discussion, see also Crawford [Q. Other important techniques for 
obtaining information on the dynamics in a vicinity of bifurcation points are perturbation methods 
which have been reviewed in detail by Iooss & Joseph [Q and Hassard et al. [fil|. These approaches 
generally provide series expansions of the non-trivial steady or periodic solutions that emerge at 
the bifurcation points. 

The paper is organized as follows. In Section ||, we develop the amplitude description for flows 
on the center manifold. Section || extends the results to describe the unfolding of the derived 
amplitude equation. The corresponding theory for iterated maps is very similar to that for flows 
and is given in Section f|. Section || illustrates the use of the general method by deriving amplitude 
equations and transformation equations for the most common simple bifurcations. It presents the 
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results as convenient tables that can be used immediately for any specific physical problem covered 
by the tables. We illustrate the use of the tables in appendix [A| by deriving amplitude equations 
and transformations for some well-known systems. 

Amplitudes as initially defined in this paper (Sections || and |]) have simple geometric signif- 
icance, common to all bifurcations for flows and iterated maps. However, for oscillatory modes, 
it is advantageous to introduce the concept of "proper amplitudes", which are complex variables 
for which both modulus and argument vary slowly. These proper amplitudes agree with accepted 
usage. It is then possible to scale time and amplitudes for all bifurcations so that the scaled am- 
plitude equations become independent of the distance from the bifurcation point to lowest order 
— we define proper amplitudes to incorporate such scaling. Scalings and proper amplitudes are 
discussed in Section ||. 

The physical interpretation of unfoldings of amplitude equations is discussed in Section |t] aiming 
at a more intuitive understanding of their meaning. Section || discusses some features of the solutions 
to amplitude equations, transformed back to the original physical space, i.e. solutions in physical 
terms to physical problems. Such back transformation is an important part of the solution and it 
follows naturally from the present method. In the final Section ^, we summarize and discuss the 
results of this paper. 



2 Amplitude equations for flows 

We start by considering a flow on R" 

x = F(x) = J-x + f(x), (1) 

with a stationary solution x = 0. Here, the vector field on the right-hand side has been split into 
linear and nonlinear parts Jx and f(x). We assume that J has r distinct right eigenvectors Uj 
and left eigenvectors u* corresponding to r critical eigenvalues Xj with zero real part (j = 1, . . . , r), 
chosen to satisfy the biorthonormality conditions 



u*-Uj = Sij, for i,j = 1, ...,r, (2) 

where <5y is the usual Kronecker delta. The r eigenvectors span the r-dimensional center subspace 
E c . As discussed in the center manifold W c can be locally parameterized by points z e P 
through a smooth transformation h : E c — > W l : 

W c = {x | x = z + h(z)}. (3) 

The geometrical interpretation of this parameterization of the center manifold is illustrated in Fig. [j]. 

The flow in W c defined by (|l|) is determined through (||) by solutions to the differential equation 
in z e E c 

(i + ^) -z = J-z + J-h(z) + f (z + h(z)), (4) 
in which I is the unit tensor. Equation (Q) may be transformed to the differential equation 

z= (l + ^)~-(j-z + J-h(z)+f(z + h(z))) (5) 
which we formally express as 

z = J-z + g(z), (6) 
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with the nonlinear part g(z) still to be determined. 

We now represent the vector z G E c by coordinates Z\, . . . , z r in the basis of critical eigenvectors 

r 

Z = J2 Z 1 U 3- ( 7 ) 

3=1 

We can then Taylor expand the relevant functions defined on E c in the coordinates Zj as follows 

f(z + h(z)) =^f P z p , (8a) 
p 

h(z) = ^h pZ P, (8b) 
p 

g(z) = ^ gpZ P. (8c) 
p 

Here 

p= (p!,...,fv) GN5, (9a) 

r 

Z P = n#, (9b) 

3=1 

and |p| = J2]=iPj — ^ for all terms of the nonlinear functions defined in (||). The coefficients 
f p , hp, and g p are vectors in M™ (C™ if any of the critical eigenvalues are complex). In terms of 
the particular coordinates Zj, the set of differential equations from (^) for the flow on the center 
manifold W c becomes 

:, A ; : ; • ]T(u*-gp)z p , (10) 
p 

which we shall refer to as the amplitude equation for ([!]) on the center manifold when the transfor- 
mation h(z) is chosen so that ( |10|) contains as few non-vanishing coefficients u*g p as possible. As 
described below, this procedure also defines the transformation h(z). 

The transformation h(z) is determined through its coefficients h p in ([sb] ) which must be obtained 
iteratively together with the coefficients f p and g p . This is possible because the center subspace E c 
is tangent to the center manifold W c at x = 0. We therefore substitute (||) and (@) in (0) to get 



(i + E E h P zP fS*) ■(£ + E gp zP 

P i=i 1 k=i 

r 

z^ku k + jJ2 h P zP + E f p zP 

k=l 



Pi=l ^ ,=1 P (n) 



The equation for the linear terms has already been anticipated in (g). For the nonlinear terms, pi] ) 
must be satisfied for each term indicated by p separately. We shall introduce a notation for the 
coefficient to z p from the product of the two sums over p on the left taken together with the last 
term on the right-hand side of dlj), namely 



*p = f P E h P' EP3 U i'S(p-p'+^)- ( 12 ) 

P' 3=1 
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Here the sum over p' is taken over all ordered sets p' = (jpi, . . . ,p' r ) for which the terms are defined, 
whereas the set Sj has Sjk as its fc'th component. 
From (fll]), we then get the coefficient g p as 

r 

gp = ( J -E^ A ^)- h p + *p- ( 13 ) 

3=1 

Thus, the coefficient g p can be eliminated if the coefficient h p of the transformation h(z) can be 
chosen in such a way that the right-hand side of (O) is zero, that is if 



-5>A-l).h p = -* p , (14) 

3=1 

is solvable for h p . By biorthogonality, ( |l4]) splits into two independent equations for components in 
E c and components in the generalized eigenspace complement to E c . In the complement subspace, 
( |l4| ) can always be satisfied. For components in E c , (14) can be solved for any component u*-h p 
of h p in the eigenvector basis u, for which Y) ■ 1 PjXj ^ Aj. For each of these components, we may 
eliminate the corresponding term u*g p of ( |lfj|) by choosing u*h p as the solution to (fl4|), namely 

< h P = T v*"''' \ ■ ( 15 ) 

In contrast, a component U*g p of g p cannot be eliminated by any choice of transformation 
h(z) if the resonance condition 

r 
3 = 1 

applies. The condition ( ]l6| ) is identical to the one appearing in a somewhat different setting in the 
traditional derivation of the normal form theorem JT(|[l7]], and the essential contents are of course 
the same in the two descriptions. The resonant terms must appear in the amplitude equation ( |ic| ) 
and are determined by ( [l3| ) once the corresponding components u*-h p of the transformation h(z) 
have been chosen. We set each of them equal to zero so any non-vanishing (resonant) coefficient 
becomes (from (|l3|)) 

u*.g p = u*.<e p . (17) 
The amplitude equation can now be written more explicitly as 

z i = A i Z|+5^«-#p>P, i=l,...,p, (18) 
p 

in which the sum for the i'th amplitude Zi is taken over all sets p that satisfy the resonance condition 
©• 

In practice, calculation of the amplitude equation ( |l8| ) does not require a decomposition of ( |l4| ) 
as we did above. For any given p, we may solve ( |l4| ) directly for all non-resonant components of h p 
in terms of the part of 3? p biorthogonal to the null space of (J — Y^j=i Pj'Ajl), so h p is completely 
determined by the set of equations 

r 

J-&Al)- h P = -Q P -*P. (19a) 

3=1 

R p h p = 0. (19b) 
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Here the projections R p and Q p are defined for any x £ K™ as 

R p -x= (I-Q p ).x = ^«.x)u„ (20) 



where the sum is taken over all integers i for which the resonance condition ( p.6|) is satisfied at 
given order p. In practice, the solution of ( |l6| ) is conveniently obtained by means of singular value 
decomposition p8| . 

We shall refer to the transformation x = z + h(z) from E c to W c with h(z) determined by 
(8b) and ( |l9| ) as the amplitude transformation associated with the (|l|). 

The problem of finding the non-vanishing coefficients of the amplitude equation ( jig] ) and the 
coefficients h p of the amplitude transformation leading to it, is solved order by order for increasing 
|p| = Y^,j=iPj- For any given p with |p| > 2, the non-vanishing components of the vector coefficient 
g p of the amplitude equation are determined by the resonant part of 3? p whereas h p is obtained 
as the solution to ( |l9| ) in terms of the non-resonant part of €> p . We note that, for any p,i, either 
u*g p or u*-h p must vanish. 

The iterative solution is possible because 3? p is determined by lower order coefficients of g p 
and h p . This fact is obvious since the sum over p' in ( fl2| ) only involves coefficients of orders 
within the range 2 < p' < |p| — 1, and since the Taylor expansion of f (x) contains no linear terms. 



Consequently, the first term, f p , on the right-hand side of (12) only contains components h p ' of 
orders 2, . . . , |p| — 1 through (j8a[ ). 

The amplitude equation fll8| ) and the linear system of equations (|is| ) together with expression 
(|l2l ) for <frp, constitute the main result of this section. Their practical use for a number of specific 
bifurcations will be demonstrated in Section [g. First we need to extend the theory to systems 
defined at a finite distance from the bifurcation point. This we do next. 



3 Unfoldings of amplitude equations 



For physical systems, one is often interested in studying the behavior of a dynamical system in a 
vicinity of a certain bifurcation point. We therefore consider a dynamical system 



x = F(x, (j,) = J-x + f (x, fi) : 



(21) 



depending on a set of parameters // € I s , in such a way that it contains the system (|l]) considered 
previously as the special case fi = 0. In physical problems, the motion in a "slow manifold" at 
[i ^ smoothly developed from the center manifold at fi = 0, is of special importance. There 
are several ways of approaching this problem. One may treat the motion in an unstable manifold 
in the same way as described in the (previous) Section 0for the center manifold with but a minor 
modification. We shall comment on this idea in Section^ 

Here we use an approach described in |H| and extend the differential equation (Ell) as 



x = J-x + f(x, /Li), 

A = o, 



(22) 



in which (x, /Li) € K n+S . Note that we treat the term J-x linear in x explicitly in ( |22| ) whereas all 
other terms are contained in f(x, /li). 

We are interested in the motion in the center manifold of the extended system (p^), which 
represents the slow motion for finite /x near fi = 0, desired for an efficient description of the 
physical problem. Unfortunately, we cannot use the results of Section || directly because the linear 
part of the vector field in (p3) need not have r + s independent critical eigenvectors. Nevertheless, 
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the structure of (22) is so simple that we can use essentially the same approach as in Section g. 
We sketch the arguments and introduce the necessary notation in order to analyze the extended 
system. 

We consider only a dependence of the vector field on the parameter fi for which (0, fi) e W l+S 
belongs to the center subspace of ( p2[ ) at (x, fi) = (0, 0). Then the center manifold of ( |22| ) may be 
parameterized by points (z, /x) £ E c x R s in the center subspace of the extended system through a 
transformation h : l n x W -> W l : 

r = {(x, / /)|x = z+h(z, fJ )}. (23) 

This means that the flow in W c defined by ( p2|) is determined by the differential equation 

* = (* + \ 3 - Z + J ' h(z ' M)) + f (Z + h(z ' M)) M) ) ' (24) 

/i, = 0. 

Clearly, this differential equation can be written in the form 

z = J-z + g(z,/i). (25) 

We express z in a basis of eigenvectors of J as z = Y^j=i z i u ji an d expand the functions f , h, and 
g in a manner similar to (j|) as 

f (z + h(z, fx), fi) = X! f Pq z P M q , (26a) 
pq 

h(z )A i) = $>p« 1 zV, ( 26b ) 
pq 

g (z,/x) = Xgpq zP ^ q ' ( 26c ) 

pq 

in which 

p= ( Pl ,..., Pr ) E%, q=(gi,...,«.)€NS, (27a) 

Note that we have introduced the short-hand notation pq for the index (p, q). Here, |p| + |q| > 1, 
with the cases |p| = 1, |q| = excluded. Since (|2^)-(|2^) each has the same structure as in section 
we arrive at the following equation for the vector g pq which must be satisfied for each separate 
order (p, q) 

r 

gpq = ( J - & V) ' h pq + *pq- ( 28 ) 
The vector <fr pq e R™ is defined by a straightforward extension of (p^), namely 

r 

* Pq = f Pq - X h p'q' XX u j -S(p-p'+*()(q-q')- ( 29 ) 
p'q' 3=1 
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Consequently, the vector g pq will vanish provided that h pq is chosen as a solution to the equation 

r 

J - ^TpjXjl) -h pq = -Qp-* pq , (30a) 

R p h pq = 0. (30b) 

Observe now that ( |30| ) has exactly the same form as jis|). We may therefore conclude that the 
unfolding of the amplitude equation in ( |l8| ) simply becomes 

£i = XiZi + ^(u*-* pq )z p /x q , i=l,...,p, (31) 
pq 

where the sum for the i'th component is taken over all sets (p, q) where resonance occurs. For 
zx = 0, only terms with |q| = contribute, and we recover the results of the previous section. For 
general /x, ( |3l] ) is an extension of ( |l8| ) which incorporates the parametric dependence on /x of the 
resonant terms. 

Equations ( p8| ) and ( psj ) thus constitute a recurrence relation for determining g pq and h pq . One 
must first determine the coefficients for |q| = for all orders p up topo + 1 where Po is the maximum 
value of |p| to be included for any q. Subsequently, the coefficients are determined for all q with 
|q| = 1 in the same way. 

On the center manifold, the functions g(z, /x) and h(z, fi) are at least quadratic in z. However, 
the unfoldings may contain nonvanishing coefficients g pq and h pq with |p| =0,1 when |q| > 1. 
In other words, the functions g(z, /x) and h(z, /x) may contain terms depending on /x which are 
constant or linear in z which now briefly discuss. For |p| = and |q| = 1, we have 

Jh 5 fc = -Qo-*o4> for k = l,...,q. (32) 

with Sk defined below (|l2|). If Ai ^ 0, the corresponding coefficient u*go4 will vanish from the 
amplitude equation provided that the coefficient u* • ho4 is chosen as 

« , U ? ^Q4 /oca 

-h 4 = — . (33) 



On the other hand, if A; = we choose u*-ho4 = by (30b), and get the following contribution 
from u* • go4 in the amplitude equation 

u*go4 = u**o4- (34) 

At the order |p| = 1, |q| = 1, we find that the following equation must hold for each ordered 
pair (<5j,4) 

:.I A,1.W Q.'l',,,,-. (35) 

For i = j, we always have u*-(J — XjTj-hs^ = 0, implying that we g et a contribution to the 
amplitude equation given by 11*^^.4 = u*-$5 3 -4 and u^h^.4 = from ( |30b| ) . In the asymmetric 
case i ^ j we have 

A^A j: u*-h aA =- U ^y , <-g 5 .4=0, (36a) 

Ai = X 3 : u* • h S]Sk =0, u* • g< 5 A = u* • $44 . (36b) 

We shall illustrate the use of these results in Section |j| and explain their meaning in Section ^. 
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4 Amplitude equations for iterated maps 

The description of the dynamics on center manifolds for iterated maps is very similar to that of 
flows presented in Section |^. Therefore we merely give a short outline of the theory. 
A discrete system in R™ is described by a map of the form 

xhF(x) = J-x + f(x). (37) 

We assume that x = is a non-hyperbolic fixed point of (|37| ) where the linear part J has r eigen- 
values Xj with = 1 (j = 1, . . . , r). The corresponding right and left eigenvectors are denoted by 
Uj and u* and are assumed to be normalized in accordance with the biorthonormal relations in (||) . 
The r-dimensional center manifold W° of the fixed point x = may be parameterized by points in 
the center subspace as 

x = z + h(z) = K(z), (38) 
where x€W c and z e E c . Inserting (||) into (|37]) yields 

z + h(z) J-z + J-h(z) +f(z + h(z)), (39) 
which shows that the discrete dynamics on the center manifold is induced by the map 

zh^ K- 1 (J-z + J-h(z) +f(z + h(z))), (40) 

in the center subspace. 

To find the map in E c explicitly, we formally express it as 

znj.z + g(z). (41) 

As in the continuous case, we express z in terms of coordinates based on the right eigenvectors Uj, 
viz. z = Y?j=i z j u j- The functions f (z + h(z)), h(z), and g(z) can then be Taylor expanded as 

f(z + h(z)) =^f P z p , (42a) 
p 

h(z) = 5> pZ P, (42b) 
p 

g(z) = ]T gp zP, (42c) 
p 

Zi i-> \iz t + ^(u*-g p )z p . (43) 



and the map ( fil| ) becomes 



When the transformation h(z) is chosen so that ( |43| ) contains as few non-vanishing coefficients 
u*-g p as possible, we shall refer to this equation as the amplitude equation for the map in (37). 

To determine h(z), we first substitute ( pl| ) for z on the left-hand side of (j39|). Since the resulting 
expression is equal to the right-hand side of (|39|), we obtain 

J z + g(z) + h(j-z + g(z)) = J z + J-h(z) + f (z + h(z)). (44) 
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Here we note that the term h(Jz + g(z)) may be written as 

r 

h(j- z + g (z)) = £ ^ liix^ + 5>;. gp ,> p ') p ' , 

P 3=1 p' 

= E h P zP fl A f +E*p zP - ( 45 ) 

P 3=1 P 

The last step defines the vector € R" which abbreviates the (rather complicated) coefficient that 

remains at order p when the term h D z p TT 7 - , A 1 * 3 is taken out. Substitution of z = y] r ,i z,u,- 

■ ■ i—i ■ 1 ~p ^ — ^ — — 

together with Q42| ) , and ( |45[ ) in ( J44| ) followed by a simple rearrangement of terms then yields the 

following expression for the coefficient vector g p 

r 

g P = (j-n^l)-hp + * P , (46) 

3=1 

where 3? p = f p \I/ p . Consequently, the coefficient vector g p can be eliminated if h p is chosen as 
a solution to the equation 

r 

(J-II^l)-hp = -*p. (47) 

3=1 

However, a component u* • g p of g p cannot be removed by any choice of h p if the resonance condition 

r 

U XP /= X ^ ( 48 ) 
3=1 

is satisfied. 

Using a series of arguments analogous to the derivation of the amplitude equation for flows given 
in section g, one easily obtains the following recurrence relation for the coefficients g p and h p of 
the amplitude equation and the associated transformation 

r 

g p =(j-n^l)-hp + *p. (49) 

3=1 



We may therefore conclude that the unfolding of the amplitude equation ([43|) becomes 

z t h-> XiZi + ^(u*-* p )z p , i = l,...,p, (50) 
p 

where the sum is taken over all resonant sets p for the i'th component. Furthermore, all resonant 
components u* • h p of h p are set equal to zero whereas non-resonant terms are removed by choosing 
the corresponding coefficient • h p of h p as 

u i 'hp = - ™^ P TpJ- (51) 

x i -ll J=1 x j 



Again, we note that the resonance condition (48) is identical to the one that appears in the tradi- 
tional normal form derivation for diffeomorphisms ||l6|,|l7|. 
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For a map depending on a set of parameters fi£l s , 

XHF(x,fi)=J'X + f(x,/j), (52) 
we use the same approach as for flows and study the extended system 

XHj. X +f(x,/l), 

(53) 

jjL i — ► fx. 

Using expansions ( p6|) and an approach completely analogous to the one discussed above, we find 
a recurrence relation 

r 

gpq= (j-II A f I )- h Pq + *P^ ( 54 ) 
3=1 

with = fp 9 — \& P9 where the definition of ^ pq is analogous to tha t of ty p . An example of 
the calculation of <P P9 for the period doubling bifurcation in Section 5J5 illustrates the use of the 



formula. In conclusion, the unfolded amplitude equation for the iterated map (|53| ) becomes 

Zi i-» \ l z l + ^(u*-* pq )z p /x q , z = l,...,r, (55) 
pq 

where the sum is taken over all sets (p, q) for which resonance occurs for the i'th component. 

The coefficients h pq of the unfolded amplitude transformation are determined as solutions to 
the linear equations 

r 

(J - JJ A^l) -h pq = -Qp-* pq , (56a) 

3 = 1 

R P h pq = 0. (56b) 

As for (|ig|), the linear system (^) is solved in practice by application of a singular value 
decomposition. 



In Section 5.5 below, we demonstrate the use of these results for the period doubling bifurcation. 



5 Amplitude equations for simple bifurcations 

In this section, we derive explicit expressions for the coefficients of the amplitude equation and the 
associated transformation h(z, jLt) to lowest non-trivial order for some simple bifurcations. For this 
purpose, we represent the vector field on the right-hand side of (|2l]) and ( |52| ) as Taylor expansions 
in some physically relevant coordinates Xj of x € R n and fij of /x € K s . In terms of x and /x the 
expansion of the vector field d2|) takes the form 



F(x,/i) = J-x + F^/x +P XM (x,/ti) + |F xx (x,x) + ^F xxx (x,x,x) + ■ • • , (57) 
in concise notation with 

F xx (x,x) - ^ a g x iXj , F xxx (x,x,x)- 2^ dx i dx i dx k XiXiXk > 

58 

V- ^F , . " " 5 2 F 

z — 1 i—l 7—1 J 
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where all derivatives are evaluated at x = and fi = 0. Here, the multilinear vector functions, 
F xx (x, x), F xxx (x, x,x), etc, are symmetric in all their arguments, so we need not specify them 
any further. We shall only consider unfoldings linear in /x, so we need not take any terms beyond 
those exhibited in ( |5^ ) into account. From part of (|57]) nonlinear in x, namely 

f(x,/x) = F^/j + F x(1 (x, /i) + iF^fx, x) + iF xxx (x,x,x) + ■ • • , (59) 

we get the coefficients f pq by first substituting the parameterization of the center manifold 

r 

x = z j u j + h pq z p /x q (60) 
i=i pq 

for x and then collecting all terms of order (p,q) (i.e. terms containing z p /x q ). In the examples 
considered here, we shall only discuss problems associated with a scalar parametric dependence [i, 
implying that we may simplify the multilinear vector functions in ( |57| ) by writing 

F xPM ,(x,x, . . . ,x,/x,/a, = F xPA1 ,(x,x, . . . ,x)fj, q . (61) 



5.1 Bifurcations at a single zero eigenvalue 

We first look at a stationary point of the flow (0) where one real eigenvalue of J is zero and assume 
that all other eigenvalues have non-zero real parts. We denote the right and left eigenvectors 
corresponding to A = by u and u* respectively. Here the center manifold and the center subspace 
are both one-dimensional, and the latter is spanned by u. There is just one amplitude z, and the 
set of integers p becomes just a simple index p. 

The resonance condition ( [l6| ) takes the form A = pX, and is satisfied for all orders p>2. So none 
of the nonlinear terms of the amplitude equation vanish in general, whereas all of the components 
u* -hp q do. At the bifurcation point, the amplitude equation ( jsil) therefore has the following form 

i = goifi + gufxz + g 2a z 2 + g 3Q z 3 H , (62) 

in which g pq = u*-g p<3 . For the lowest order contributions to the amplitude equation and the 
transformation h(z, ^), we need first calculate 3> pg from (|2g|). To obtain the coefficients f P9 , we 
substitute 



x = uz + h i/i + hn^z + h 2 o2 2 + h 30 z 3 + • • • (63) 



in the expansion (pj 



f(x) = Ff,fi+ (F XM -u + F xx (u, h i))/x2 



iF xx (u, u)z 2 + (F xx (u, h 20 ) + |F xxx (u, u, u))* 3 



(64) 



Thus, 



foi = P M , (65a) 

fix =F XM -u + F xx (u,h i), (65b) 

f 2 o = iF xx (u,u), (65c) 

f 3 o = F xx (u, h 20 ) + £F xxx (u, u, u). (65d) 

Second order in z (p = 2 and q = 0) 

At second order, there is no contribution from the sum over p' and q' in (p9|), implying that 



*2o = f 20 = 2-F xx (u,u). (66) 
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We therefore get the second order coefficient of the amplitude equation from ( |3l| ) as 

.920 = u*-* 20 = ±u*-F xx (u,u), (67) 

and a linear equation for the coefficient vector I120 from ( ^o|) 

Jh 20 = -±Q-F xx (u,u), (68a) 
u*-h 20 = 0, (68b) 

with the projection Q defined through Qx = x (u*-x)u. 
Third order in z (p = 3 and q = 0) 

Once g2o and h 2 o have been obtained from (|67| ) and (|6|) respectively, we may derive the third 
order contributions using (H) for 3>3o 



*30 — f30 — 2.920 n 20 j 

= F xx (u, h 20 ) + ^F xxx (u, u, u) - 2 5 2oh2o- (69) 

We observe that u*-h 2 o = 0, and conclude from the above expression that 

.930 = u*-* 30 = u*-F xx (u,ti 20 ) + ±u*F xxx (u,u,u), (70) 

whereas the coefficient vector h 30 is determined by the linear equation 

Jh 30 = Q(F xx (u, h 20 ) + ^F xxx (u, u, u))+2 52 oh 2 o, (71a) 
u*-h 3 o = 0. (71b) 

Following this procedure, the iterative solution for g p o and h p o can be developed to still higher 
orders. 

Unfolding terms (p = and 9=1) 

Here 3>oi = foi- From (|65|), we get the first contribution to the unfolding as 

,9oi=u*-F M . (72) 

In order to find h i, we must solve the linear equations 

J-hoi = Q F M , (73a) 
u*-h O i=0. (73b) 

Unfolding terms (p = 1 and 9=1) 

From (p9|), we find that $n = fn — 2901I120. Having obtained both goi and hoi, we may finally 
find the contribution gn to the unfolding of the amplitude equation ([62]) as 



.911 = u*-*n = u*-F XjU -u + u*-F xx (u,h i), (74) 

since u*-h 2 o = 0. The coefficient vector hn can then be obtained by solving the linear system of 
equations 

J-hn = -Q(-F XM -u + -F xx (u,h i))-2.9oih 2 o, (75a) 
u*-hn = 0. (75b) 

We now specialize the results for the amplitude equation to three important bifurcations at a 
single zero eigenvalue, including only the terms at lowest order in z for q = 0, 1. For the terms 
linear in /i, we introduce the simpler notation a p = g p i, whereas we shall use g p = g p o for terms 
independent of [i. 
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5.1.1 Saddle-node bifurcation 

If (To 7^ and g 2 7^ 0, a saddle-node bifurcation occurs at /i = 0, and the corresponding amplitude 
equations takes the form 

i = (7o/i + g2Z 2 , (76) 



with 



.92 = ±u*-F xx (u,u), (77a) 
<7 = u*-F M . (77b) 



5.1.2 Transcritical bifurcation 



If (To =0 but <ti and 172 7^ 0, a transcritical bifurcation occurs at n = 0, for which the amplitude 
equation is 

z = <Ji/.iz + g 2 z 2 , (78) 



with gi given by (77a) and 

01 =u*-F X(U -u + u*-P xx (u,h i). (79) 

5.1.3 Pitchfork bifurcation 

Finally, if <tq = and .92 = 0, but o\ 7^ and 33 7^ 0, a pitchfork bifurcation is realized at /i = 0. 
Here the amplitude equation becomes 

i = aifiz + g3z 3 , (80) 

where o\ is given by (|7|) and 

53 = u* F xx (u, h 20 ) + ±u* F xxx (u, u, u). (81) 

The results derived for these three bifurcations are summarized in Table [j] which also exhibits the 
equations for the coefficients of the amplitude transformation. 

5.2 Bifurcations at A = 1 for iterated maps 

The results derived above for flows require little modification in order to give the coefficients of 
the equivalent amplitude equations for saddle-node, transcritical, and pitchfork bifurcations in the 
iterated map ( |37| ) (at a fixed point where one eigenvalue of J satisfies A = 1 and all other eigenvalues 
have A 7^ 1). In fact, the only necessary modification is to 

replace J by J — I (82) 

everywhere in (^), (f7l|), (|73|), and (|75|). For this reason, we do not show the details of the calcula- 
tions for the iterated map here, but simply give the results for the unfolded amplitude equations: 

z I— > (To/- 1 + z + 52^ 2 (saddle-node), (83a) 

z I— > (1 + a\{i)z + g 2 z 2 (transcritical), (83b) 
z 1— > (1 + (Tin)z + g3Z 3 (pitchfork) . (83c) 



Here the coefficients cjq, <ti, g 2 , and 33 are given by ( ]77b| ) , ([79]), ( 77a ), and ( |8l] ) with the linear 
equations ( |68| ) , (^l]) , ( f73] ) , and ( |75| ) modified in accordance with (|82 ) . The results derived for 
the saddle-node, the transcritical, and pitchfork bifurcations for iterated maps are summarized in 
Table § 
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5.3 Hopf bifurcation 

We now consider a stationary point of the flow ( ^l|) at which the linearization J has two complex 
conjugate, pure imaginary eigenvalues Ai = A2 = iu>o, and all other eigenvalues have non-zero real 
parts. We denote the right and left eigenvectors of J corresponding to Ai by ui and u*, and those 
of A2 by U2 and u 2 , as in the general theory — and choose U2 = Ui normalized according to (||) . 
Consequently, we may simplify the notation by defining u = Ui and u* = u*. By this choice, the 
two-dimensional center subspace is spanned by any linear combinations of Re u and Im u, whereas 
the motion in the two-dimensional center manifold is determined by an amplitude equation which, 
as a straightforward application of the general theory, is formulated in terms of the two complex 
eigenvalues and complex conjugate amplitudes z,z with z = zu + zu G E c . 

The ordered set of integers p now has two components, p = (pi,p2). The resonance condition 
for any p is satisfied for either of the two components if 

P1A1 +p 2 A 2 = (pi -p2)Ai = Ai (84a) 



or 



P1A1 +P2A2 = (p2 -pi)^2 = A2. (84b) 

Thus, the resonance condition is fulfilled for the first component if pi = pi + 1, and for the second 
component if P2 = Pi + 1 provided that p\ + P2 > 2. The amplitude equation therefore takes the 
form 

z = iu> z + g W iLiz + g 2W \z\ 2 z + g 32 o\z\ 4 z -\ , (85) 

where g ljk = u*-g ijk . 

To derive expressions for the coefficients of (|S5|), we first determine the constant vectors f pq by 
substitution of 

x = uz + uz + h 2 oo-z 2 + hno|z| 2 + h 2oz 2 + h oi^ + h W i^z + honfiz + ■ ■ ■ (86) 

in (|59|). Collecting separate orders in z, z, and /i, we find the following contributions for terms of 
zero order in [i 

f 2 oo = ±F xx (u,u), (87a) 
fiio = F«(u,u), (87b) 
f300 = F xx (u, h 20 o) + gF xxx (u, u, u), (87c) 
f 2 io = F xx (u,hno) + F xx (u,h 20 o) + iF xxx (u,u,u). (87d) 



whereas the terms linear in \x are 



fooi — F M , 



fioi =F XM -u + F xx (u,h 00 i), (88b) 

Notice also that $ p o = f P o at second and third orders in |p| because the first non-vanishing 
coefficients g p o come at third order in |p|. Although all second order coefficients g p o vanish, we 
must still derive the coefficients h p o of the transformation h at second order since these are needed 
to obtain the higher order coefficients g p o as well as for determining the amplitude transformation. 
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Second order terms in z (pi + p 2 = 2 and q = 0) 
From (|87|), we find 



*2oo = *02o = ^F xx (u, u), (89a) 
*uo = Fxx(u,u). (89b) 



Since there is no resonance at second order, the right-hand side of ( |30a| ) is simply — 3? p o, and 
we get using Ai = A 2 = iujq 

(J-2iw I)-h 20 o = -5F xx (u,u), (90a) 
Jhi 10 = -F xx (u ! u), (90b) 

whereas ho2o is determined by the relation h.020 = h 2 oo. 
Third order terms in z (pi + p 2 = 3 and q = 0) 

Once ( pCf ) has been solved for the second order coefficient vectors h p o, we immediately get the third 
order terms <fr p o explicitly from ( p7| ) 

*300 = F xx (u, h 20 o) + ^F xxx (u, u, u), (91a) 
*2io = F xx (u,hno) +F xx (u,h 200 ) + ±F xxx (u,u,u), (91b) 

*120 = *210, *030 = *300- (91c) 



The third order nonlinear coefficient 5210 of the amplitude equation (85) then becomes 



ff2io = u*-*2io = u*-F xx (u,hno) + u*-F xx (u,h 20 o) + 5U* F xxx (u, u, u). (92) 

Notice that expressions like ( |9l| ) for $ p o can be written down immediately in terms of the coefficients 
h p o and g p o, but to get explicit expressions, equations like ( |90| ) must first be solved for the lower 
order coefficients. In other words — to get explicit expressions, one must work order by order 
starting with the second order. 

For the amplitude equation (|s5r) to third order, we need go no further. However, to obtain an 
explicit expression for the transformation h(z) to third order we must solve the linear equations 

(J - 3iw I)-h 3 oo = -*300, (93a) 
(J-ia; I)-h 2 io = -Q-*2io, (93b) 
u*-h 210 = 0, (93c) 

h 30 = h 30 o, hi 20 = h 2 i , (93d) 

where the projection Q is defined as 

Q-x = x (u*x)u. (94) 
Unfolding terms (pi = p 2 = and q = 1) 

First, the coefficient vector hooi, must be found by solving the linear equation 

J-hooi = -F M , (95) 
which is non-singular since Ai = A 2 = iujq 7^ 0. 
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Unfolding terms (pi = 1, P2 = and q = 1) 

We first calculate 3>ioi as 

*ioi = fioi - 2h 20 o(u* -gooi) - 2hno(u* -gooi) (96) 
= fioi =F XjU -u + F xx (u,h 00 i), (97) 

since gooi = 0. Having obtained hooi, the first non-trivial coefficient of the unfolded amplitude 
equation becomes 

gioi = u*-$ioi = u*-F XM -u + u*-F xx (u,h 00 i), (98) 

whereas the linear equation which determines the corresponding coefficient vector hi i in the trans- 
formation becomes 

(J-icJoI)-hioi = -Q-*ioi, (99a) 
u*-hioi=0, (99b) 
hon = hioi, (99c) 

We now simplify the notation using the fact that one index, p — p\ + P2 say, is sufficient to describe 
the resonant terms since p\ = (p+l)/2andp2 = (p— 1)/2 for the first component. Setting o\ =3101 
and 33 = 32101 we finally conclude that the unfolded amplitude equation becomes 

z = (icd + (7i/i)z + g 3 z\z\ 2 , (100) 

with 

g 3 = u*-F xx (u,hno) + u*-F xx (u,h 20 o) + 5U* F xxx (u, u, u), (101a) 

cri =u*-F XM -u + u*-F xx (u,h 00 i). (101b) 

The results derived for the Hopf bifurcation for flows are summarized in Table ||. 



5.4 Neimark-Sacker bifurcation 



The equivalent situation where an iterated map (j32]) exhibits a Neimark-Sacker bifurcation corre- 
sponds to a fixed point where two complex eigenvalues of J satisfy Ai 



A 2 



^ (|Ai| = |A 2 | = l). 

For simplicity, we exclude the strong resonance cases where 60 = ^ for k — 1, 2, 3, or 4. The 
unfolded amplitude equation and the associated coefficients for this problem is easily obtained by 
only a few modifications of the results derived above for the Hopf bifurcation for flows. The only 
required modification is to 



replace (J — i /cwqI) by (J 



I) 



(102) 



everywhere in ( |90| ) , ( J93| ) , ( |95[ ) , and (|99j) . We shall not show the details of the remaining calcula- 
tions here but merely present the main result for the amplitude equation for the Neimark-Sacker 
bifurcation, namely 



z ^ ( e iflo + a lf x)z + g 3 z\z\' 



(103) 



with the coefficients 33 and 01 given by (101a) and ( |l01b ) respectively with the linear equations 
(|90|), (|93|), and ( |9^ ) modified according to (102). The results for the amplitude equation for the 
Neimark-Sacker bifurcation are summarized in Table B 
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5.5 Period doubling bifurcation 

To illustrate the principles of deriving amplitude equations for iterated maps, we consider a fixed 



point solution x = of ( |53| ) where one eigenvalue of J satisfies A = — 1 at // = whereas all other 
eigenvalues are assumed to have modulus different from unity. The right and left eigenvectors are 
denoted by u and u* respectively corresponding to a one-dimensional center manifold described by 
a single real amplitude z and a simple index p. 

The resonance condition ( fl8|) i s simply A = A p , which is satisfied for all odd p > 3. Correspond- 
ingly, the amplitude equation (|50D takes the form 

z i ^ (gufjL - l)z + g 3a z 3 + g 5Q z 5 H , (104) 



with g pq = u* g 



pq- 



To determine relations for the coefficients of (104), we first insert 



x = uz + h 2 oz 2 + h 30 2: 3 + hoi// + hn/iz H (105) 

in ( |59| ) and find the following expressions for f2o, f3o, foi, and fn 

f 20 = ±F xx (u,u), (106a) 

f 3 o = F xx (u, h 20 ) + ±F xxx (u, u, u), (106b) 

fci =F M , (106c) 

fn =F XM -u + F xx (u,h i). (106d) 



Now, to derive expressions for the constant vectors <fr p o defined implicitly b y (|44| )-(|54|), we note 
that u* -g 2 o = 0, and conclude that * p o = for p < 3 since u* -g p o first enters (|45|) at fourth order. 
This therefore implies $20 = f 2 o and €>3o = f3o- 

Second order terms in z (p = 2 and q = 0) 

To find h 2 o, we have 

*20 = iF xx (u,u). (107) 
There is no resonance at second order, but we must still solve the regular linear system 

(J-I)-h 20 = -±F xx (u,u), (108) 
for h2o in order to obtain an explicit expression for $30. 

Third order terms in z (p = 3 and q = 0) 

For $30, we have 

*30 = F xx (u, h 20 ) + ^F xxx (u, u, u), (109) 
so once the linear equation for h2o has been solved, we may determine the resonant third order 



coefficient in ( [104] ) as 

330 = u*-* 30 = u*-(F xx (u,h 20 ) + |F xxx (u,u,u)). (110) 

The third order term h3o of the transformation h(z) can then be found by solving the equations 

(J + I)-h 30 = -Q-* 3 o, (111a) 
u*-h 30 = 0, (111b) 
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where the projection Q is defined in (|20|), i.e. Qx = x (u*-x)u. 
Unfolding terms (p = 0, 1 and q = 1) 



To find 3>oi and <&n, we use & pq = f pq — ^ pq and take f pq from (106). To determine any contribution 
from ^ pq , we substitute the right-hand side of ( 104) for z in the expansion of h(z, fx). Now, substract 
h(— z, fi) and order the rest according to (p,q). The coefficient of z p pfl is then <& M . We may than 
conclude for the period doubling bifurcation that $oi = *n = 0. Consequently, 

*oi = F M , (112a) 
*n =F XM -u + F xx (u,h i). (112b) 

To determine the coefficient vector hoi , we must solve the linear system 

(J-I)-hoi = -F M , (113) 

which always is regular since J by assumption has no eigenvalue A = 1. Having obtained hoi, we 
find the non-trivial contribution to the unfolding of the amplitude equation as 

u*gn = u*F x/i u + u*F xx (u,h i), (114) 

and the coefficient vector hn by solving the linear system of equations 

(J+I)-h u = -Q-(F x/1 -u + F xx (u,h i)), (115a) 
u*-h u =0. (115b) 

In conclusion, we arrive at the following unfolded amplitude equation for the period doubling 
bifurcation 

(aui-l)z + g 3 z 3 , (116) 

where 

.93 = u*-F xx (u,h 20 ) + ±u*-F xxx (u,u,u), (117) 
«n = u*-F XM -u + u*-F xx (u,h i). (118) 

The results for the period doubling bifurcation are summarized in Table || 



6 Scalings and proper amplitude equations 

The amplitude of a harmonic oscillation is half the difference between the maximum and minimum 
values of the oscillating function. This basic concept of an amplitude can be generalized to a set 
of (possibly complex) Fourier amplitudes for anharmonic periodic oscillations. For the aperiodic 
(transient) oscillations that can occur n ear a Hopf bifurcation, we can introduce a slowly varying, 
complex amplitude w related to z of ( |100| ) by extracting the time dependence associated with 
harmonic oscillations with (vanishingly) small amplitude at the bifurcation point. Specifically, we 
can extract a phase factor e luJat from z(t) to get a time dependent amplitude modulating the "basic" 
oscillations described by the phase e 1 " '. 

We shall introduce the notion of a "proper" amplitude to account for this and similar slow mod- 
ulations. At the same time, it is convenient to account also for the important scalings characteristic 
of a Hopf bifurcation. Thus, the amplitude of the limit cycle oscillations grows as y/JI for small fj,, 
and the characteristic times in a vicinity of the stationary point as well as near the limit cycle are 
both proportional to /i close enough to the bifurcation point. 
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We therefore introduce a proper (scaled) amplitude w(t) as a function of a scaled time r by the 
definition 



w(t) = -^—exp(-iw t)z(t), 
r = [it. 



(119) 
(120) 



Substitution of these two expressions in (10C) then leads to the proper amplitude equation for w{t) 



dw 
d7 



o-iw +g 3 w\w\' 



(121) 



which traditionally is called the Stuart-Landau equation |l^,^o|. It has the important property 
of being independent of the distance from the bifurcation point (which of course only applies at 
the lowest non-trivial orders in z and [i of the original amplitude equation) . So the Stuart-Landau 
equation is universal for any flow exhibiting a Hopf bifurcation — characterized by just two pa- 
rameters <j\ and §3. Note that the scalings leading tojl21) are dictated by the form of ( |100| ) . The 
Stuart-Landau equation (121) can be derived from ( lOOD w ithout any knowledge of the properties 
of the Hopf bifurcation othe r tha n those contained in (100). 

Our remarks leading to (121) are meant solely as motivation of the idea of scaling and are not 
needed for the derivation itself. On the contrary, we can infer the scaling properties for the limit 
cycle and the critical slowing down from the amplitude equation. This feature distinguishes the 
present derivation of (proper) amplitude equations from those based on multiple time scales; in the 
latter case the scalings are basic assumptions of the actual derivation of the amplitude equation 
whereas these follow naturally from the derivations presented here. 

For non-oscillatory modes, the amplitude equations at lowest nonlinear order can also be scaled 
to yield a proper amplitude equation which also will be universal for any flow exhibiting the bifurca- 
tion in question. In fact, similar proper amplitude equations are easily derived for the saddle-node, 
the transcritical, and the pitchfork bifurcation. For the equivalent amplitude equations of an iter- 
ated map, one can show that the scaling t = [it is replaced by an explicit factor \i in the scaled 
amplitude equation. 



7 Interpretation of the unfoldings 

The significance of the description of the motion in the center manifold in terms of an amplitude 
equation is obvious from the geometrical interpretation illustrated in Fig. 0. But the meaning of 
the unfolding is perhaps not so clear. 

To understand terms like the coefficient vector h pq in (|30| ) for |q| ^ 0, we look first at the case 
|p| = 0, |q| = 1 corresponding to ([32|)-(|34]), These equations are associated with the change of the 
stationary point x s (/x) with the parameter fi as described by the equation 

F(x.(/i),/i) =0. (122) 

The change from (x s , fi) = (0, 0) is to first order in fi determined by the equation 

-•^(J.— + F,).^0 (123) 

in which all de rivat ives are evaluated at (x, fi) = (0,0). For finite fi, the expression inside the 



parentheses in (122) must vanish, and we find 

<9x 
<9/x 



J-^ = -F M , (124) 
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which corresponds to (32) for the fc'th component of fi. 

Thus the component of hos k in the range of J determines the change in 
point x 8 (/n) as 



of the stationary 



*0) = J2 h os h 



(125) 



fc=i 



A change with fi of the motion in the nullspace of J is determined by the contribution to the 
amplitude equation ([n]) provided the derivative F M has a non-vanishing component in the null- 
space of J. In contrast, the transformation h will have no component in the nullspace. 

We see the effects of the |p| = 0, |q| = 1 contributions to the unfoldings for the single zero 
eigenvalue in Section 5A and for the Hopf bifurcation in Section EL3. Thus, the change of the 
stationary point is given to first order in fi by ( f73| ) and (^5|). In this case, there is only a contribution 
to the amplitude equation for the saddle-node bifurcation where F M has a non- vanishing component 
in the nullspace of J. Here the non- hyperbolic stationary point splits in two which cannot be 
accounted for as a smooth displacement of the stationary point under the unfolding — a term must 
therefore appear in the amplitude equation. For the transcritical, pitchfork, and Hopf bifurcations, 
on the other hand, the change of the stationary point with [i can simply be accounted for by the 
smooth transformation h(z, /x). 

At the order |p| = 1, |q| = 1 the unfolding can be understood in terms of perturbation theory 
for the eigenvalue problem. Consider the change of the Jacobian matrix at the stationary point x s 
to first order in fj, 



3(H) = 3 + FjcO.m) + F xx (-,x s ), 



(126) 



in which x s (to first order in fi) is given by (125). For short we write (126) as 

3{n)=3 + S3 (127) 

and consider the eigenvalue problem 

(J + 63) ■ (uj + Su 3 ) = (\j + SXj)(uj + 5 Uj ). (128) 

To lowest order (linear in fi), we find from perturbation theory J2l]] that the change SXj of the 
eigenvalue Xj is 

SXj =u*-63-Uj, (129) 
whereas the change Suj of the eigenvector Uj is determined by the linear equation 



(J — XjT) ■ Suj = — (I — UjU*j ) ■ S3- Uj , 



u* -Suj = 0. 



From (126), we find for 63 Uj 



S3 Uj = F XM (uj,/x) + F xx (uj,x 8 ). 



(130) 



(131) 



By comparison of (|l30| ) and (|13l] ) with ( |35|) , we identify hg^ as the contribution to the change 
of the j'th eigenvector per unit change of [ik- Similarly, we identify the coefficient u*-g^,5 fc in the 
unfolded amplitude equation as the derivative of the j'th eigenvalue with respect to 
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The effects of the |p| = 1, |q| = 1 contributions can be inferred directly from the specific 



bifurcations studied in Sections 5.1 and 5.3. For a scalar parameter /i, we have <&SjS k = *&5ji 



implying that (131) reduces to 

53-Uj = F^-uj +F xx (u 3 -,h i). (132) 



The right-hand side of this expression has the same form as ( |65bD and ( |88b| ) for a single zero 
eigenvalue and for a Hopf bifurcation respectively. For a scalar parameter, we also note that the 
change o f th e st ation ary point is simply hoi/i corresponding to first order terms hgi/x and hgoi/x in 



Sections |5.l| and |5.3| respectively. 

We can summarize these results for one critical component with eigenvalue A, right eigenvector 
u and amplitude z in simplified notation (disregarding other critical modes) by grouping terms as 
follows 

i = goifi + (A + g\\\i)z + • • • , (133a) 
x = h i/i + (u + h u /i)z + • • ■ , (133b) 

where '• • • ' indicates linear contributions from other critical components plus nonlinear terms. 



Thus, the result (133) partly accounts for the change of the stationary point to and changed 
eigenvalue A + gnfj, and eigenvector u + hn/z evaluated at the new stationary point. It seems 
reasonably to assume that other terms can be interpreted as corrections changing the amplitude, 
but we shall not pursue the interpretation any further. 

The interpretation indicated above suggest that for finite fJ, one may obtain an efficient de- 
scription of the dynamics by using the formulas obtained for the center manifold using the (nearly 
critical) eigenvalues and eigenvectors as well as all coefficients calculated from the vector field at 
x s (/x) for the finite value of /x. 

In the important special case where no eigenvalues have positive real parts at the bifurcation and 
all the critical eigenvalues and no others acquire positive real parts at the finite /x, we can appeal 
to the stable manifold theorem and use essentially the same line of reasoning as in Section || for the 
description of the motion in the unstable manifold. We may still classify the terms to include in the 
amplitude equation and those to include in the associated transformation according to whether or 
not the resonance condition is satisfied at the bifurcation point. Although it is in general possible to 
eliminate all nonlinear terms from the amplitude equation, one should not do so, but instead include 
precisely those that appear at the bifurcation. Then the amplitude equation and transformation 
change smoothly with /i., and the formulas derived here for the center manifold can be applied to 
the unstable manifold with all quantities calculated at x s (/x) for finite fi. 



8 Physical relevance of the amplitude equation 

An amplitude equation describes the kind of dynamical behavior (e.g. chaos) that can be observed 
in any system to which it applies. It can therefore be studied in its own right as a representative of 
a large class of dynamical systems. 

For a specific system, the amplitude equation has its coefficients fixed by the physical system. 
Some or all of these coefficients can often be determined from a summary knowledge of the system, 
such as characteristic (invariant) exponents. Furthermore, since the actual (physical) system is 
topologically equivalent to that determined by the amplitude equation, the number of solutions 
that branch off at the bifurcation point and the stability of the bifurcating solutions (the super- and 
subcritical nature of the bifurcation) are additional examples of characteristics that are determined 
directly by the properties of the amplitude equation. 
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However, to describe a real physical system in more detail in terms of amplitudes, the amplitude 
equation is not sufficient. The transformation ( |60| ) between the description in terms of amplitudes 
and the actual physical variables is essential. Fortunately, that transformation is obtained together 
with the amplitude equation with little extra work, since the coefficients h pq of the transformation 
must be calculated in order to determine the amplitude equation anyway. 

For a given amplitude equation, we may therefore study the dynamics by either analytic or 
numerical methods and then use (pQ ) to transform the dynamics back to the physical space of 
the original problem. For all the amplitude equations derived in Section |[ an analytic approach 
is possible. The calculation of the limit sets of these amplitude equations is straightforward. The 
results and their relation to the physical space is summarized in Tables and Table f| respectively. 

The advantage of using an amplitude approach to a physical problem is that the essential 
dynamics usually takes place in a subspace (the center subspace of amplitudes) of much lower 
dimension than the full dimension of the dynamical system. So whether one has an experiment 
or a model, it is often advantageous to solve or analyze the problem in terms of amplitudes. For 
such work, a quantitative approach like that described in the present paper is necessary. As an 
example we mention a quantitative experimental analysis of oscillato ry c hemical reactions near a 
supercritical Hopf bifurcation relying on the Stuart-Landau equation (121) and the basic geometry 
of the amplitude transformation |22|, and a study of chemical waves in a complex chemical 
reaction-diffusion system 0|. 



9 Summary and discussion 

The main result of this paper is a general and systematic method of deriving amplitude equations 
(or normal forms) for flows and iterated maps near a local bifurcation. In addition, the method also 
provides the transformations that connect the amplitude description with the "real world" described 
in terms of appropriate physical variables, the amplitude transformation. In particular, the method 
provides an explicit recurrence relation (|||) with ( p9| ) for flows, and (|J) for iterated maps in the 
case of semisimple critical eigenvalues. 

The method is completely transparent: the amplitudes representing the instantaneous state 
of a system and the amplitude transformation to a physical description of the state have simple 
geometric interpretations. The use of the method for a given type of bifurcation is straightforward 
and unambiguous, and can in principle be carried out to any desired order of nonlinearity. 

The method is efficient in terms of necessary computational efforts as will become clear from a 
summary of the steps of the method which we now present: We first have to find the eigenvalues A^ 
and the associated eigenvectors and u* satisfying the condition RcA^ = (for flows) or |Ai| = 1 
(for iterated maps). From this, we immediately get the structure of the amplitude equation and the 
amplitude transformation h(z,/x) through a classification of terms u*g pq and u*h pq depending 
on whether or not the resonance condition 

r 

^PjAj = Xi (for flows), (134a) 

r 

\^ = A, (for maps) (134b) 

j=i 

is satisfied for the i'th component (amplitude Zi) of the term indicated by the order (p,q). The 
amplitude equation includes all resonant terms and no others, whereas the amplitude transformation 
contains non-resonant terms only. 

The next step is to calculate the constant vectors <fr pq in terms of the Taylor coefficients f pq , 
g pq , and h pq . This step is straightforward. For simple bifurcations (as considered in this paper), 
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the necessary calculations are easily done by hand, whereas more complex problems of higher 
codimension is a simple task for software tools such as Mathematica and Maple p4),|25||. 

The coefficients g pq (for the amplitude equation) and h pq (for the amplitude transformation) 
are now obtained order by order in an iterative procedure. At any order (p,q), any non- vanishing 
component u* • g pq in the amplitude equation is given explicitly by 

<gpq = <*pq- (135) 

This term only appears if the resonance condition is fulfilled. To obtain the coefficient h pq , one has 
to solve the system of linear equations 

r 

J - y^pjXjl) - hpq = -Q p -<P pq (for flows), (136a) 

i=i 

r 

J ]J x 7 T ) ' h pq = -Qp-*pq ( for ma P s ) ( 136b ) 

determining the non-resonant components of h pq , with the auxiliary condition 

Rp-hp q = 0, (137) 



which ens ures that the resonant components of h pq vanish. Equation (|136a|) or ( |136b| ) together 



with (137) determine h pq completely. 

This is all that is needed to formulate a given physical problem in terms of an amplitude 
equation. Notice carefully that the above procedure does not require that the Jacobian matrix 
has been transformed to block form — a feature that greatly simplifies the computations. The 
only non-trivial parts of the computations are the determination of the right and left eigenvectors 
associated with the bifurcating eigenvalues and the solution of the systems of linear equations ( |136| ) . 

Another significant result of the present paper is the compilation of general results for several 
types of bifurcations as summarized in Tables Q-^. This part of the work serves two important 
purposes. First we use it to demonstrate how the general procedure (sketched above) works in 
practice for flows and iterated maps. Secondly, it provides explicit expressions that can be used 
immediately for any specific problem. Such actual use of the tables is illustrated in Appendix |a] 
where the calculation of amplitude equations and associated transformations is demonstrated for a 
few specific problems, using the formulas of the tables. 

As can be seen from the specific bifurcations discussed in Section H and the examples of the 
Appendix, the machinery is quite easy to use despite the array of expansions involved. In prac- 
tice, these expansions are finite or can be truncated so that only few terms need be considered. 
Usually, only a few critical modes are involved. Thus the amplitude equation has few components, 
and in general one needs only include its first few terms to get a satisfactory description of the 
dynamics. The associated amplitude transformation can usually be truncated after the highest or- 
der included in the amplitude equation. We also emphasize that considerable simplification occurs 
in many applications as a result of the form of the vector field or map. For example, this is the 
case in chemical kinetics based on elementary reactions, which are at most of the second order. 
Consequently, the Taylor expansion $5l\ ) terminates at low order because higher order derivatives 
vanish. As mentioned previously, if the derivations by hand become too laborious, one can always 
carry it through by resorting to software package like Mathematica or Maple which are capable of 
performing symbolic computations 

The main application of the results of the present paper is for determining the time evolution 
of dynamical systems near a local bifurcation. However, the derivation of an amplitude equation 
also provides an explicit expression for the center manifold in the simplest possible way. 
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In addition, the formulas of Tables |]^4| provide a general tool for determining numerical esti- 
mates of branchings of stationary, fixed points, and periodic solutions of flows and maps when a 
bifurcation is traversed. Such estimates are typically needed in continuation problems as imple- 
mented in the numerical software packages CONT pp. 302-361] and AUTO Q. Estimates 
based on the solutions of amplitude equations are more accurate than those traditionally used in 
terms of tangents of the intersecting branches at the bifurcation point |^9| . These do not include 
the curvature of the branches at the bifurcation and equally important — they do not determine the 
stability character (as super- or subcritical) of the bifurcation point. All of these issues are generally 
important when calculating delicate solution and bifurcation diagrams. For details re gar ding the 
numerical location of bifurcation points of the type discussed in this paper, we refer to |2g-|3l| and 
references therein. 

Notice also that the method presented here is generally valid also for stable and unstable man- 
ifolds. Except for accidental resonance, the differential equation describing the dynamics on these 
manifolds can in general be linearized to any order implying that all terms are non-resonant and 
therefore contained in the associated transformation. 

As already mentioned in Sections || and |J, we emphasize that the theory for amplitude equations 
can be used directly to perform reductions of flows and iterated maps to normal form. Traditional 
approaches, where the normal form theorem is applied directly to a given problem, involves the 
determination of Lie brackets and solutions of the associated homological equation. As noted 
by Arrowsmith & Place fll^ , such calculations are often complicated and rather lengthy. The 
approach considered here, however, determines the normal form and the associated (normal form) 
transformation simply by solving systems of linear equations. From practical experience, we find the 
approach considered here much more convenient compared to traditional normal form reductions. 

In this paper, we have worked out explicit formulas for a number of simple bifurcations. The 
derivation of these serve as an illustration of the general method presented but also provides a 
collection of results than can be useful for applications to specific systems: it may indeed save much 
time. 

However, the method of deriving amplitude equations and associated transformations is abso- 
lutely straightforward also for more complicated local bifurcations (of codimension-two and higher) . 
We have derived explicit amplitude equations and transformations for some of these higher bifur- 
cations for flows and maps with the method described here. In this connection, it is natural to 
treat also bifurcations where the geometric multiplicity of the critical eigenvalues is smaller than 
the algebraic multiplicity. Such bifurcations can be treated in terms of generalized eigenvectors by 
straightforward modification of the present method. However, the general case requires a compli- 
cated notation, so we have not made the extension to generalized center eigenspace in the present 
paper. At the time of writing, we have worked out explicit expressions for coefficients of ampli- 
tude equations and transformations for a Bautin bifurcation (where a Hopf bifurcation changes 
from super- to subcritical and a quintic term is needed), for a Bogdanov-Takens bifurcation, for 
a fold-Hopf bifurcation, and for a "Hopf-Hopf" bifurcation (where two pairs of complex conjugate 
eigenvalues pass the imaginary axis) with and without strong resonances. For iterated maps, we 
have worked out the four cases of strong resonance for the Neimark-Sacker bifurcation. 

In a future paper, we shall show how an amplitude equation for a fold-Hopf bifurcation can 
be applied on a quantitative level to describe a realistic reaction-diffusion model system for the 
Belousov-Zhabotinsky reaction under circumstances where the complex Ginzburg-Landau equation 
fails entirely. In fact, this observation suggests that such codimension-two descriptions may be 
much more useful than the codimension initially suggests. In particular, that paper will illustrate 
how results for amplitude equations and their associated transformations can be used to obtain a 
description in physical terms that may be compared quantitatively with the full model system and 
experimental data. 
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A Examples 

In order to illustrate the theory discussed in the previous sections and demonstrate the use of 
the formulas of Tables Q-Q, we discuss three simple examples. They will show the fundamental 
steps involved in calculating the amplitude equation for various bifurcations in actual systems. We 
emphasize that the use of amplitude equations for the particular examples discussed here has been 
chosen to illustrate the principles of the method for some simple bifurcations — not because it 
necessarily provides the computationally shortest path to the solution. 



A.l The Lorenz equations 

To illustrate the derivation of amplitude equations for bifurcations involving a single zero eigenvalue, 
we consider the Lorenz equations j3^] 

s(x 2 - Xl) 

rxi — x 2 — X1X3 ] , (138) 
-6x3 + X1X2 









-( 







and choose the term r as the bifurcation parameter. One easily sees that ( |138|) admits the trivial 
solution (0,0, 0) T with the corresponding Jacobian 



J(r) = r -1 , (139) 




which has a single real eigenvalue A = for r = 1. This implies that the Lorenz system has a 
one-dimensional center manifold at the stationary point (0,0, 0) T for r = 1. The right and left 
eigenvectors u and u* of J = J(l) are 

u=(l,l,0) T , (140) 
L -(1,8,0). (141) 



1 

Second and third order terms 

Since the coefficient oq in Table |j] satisfies oq = u*-F r = u * ■ ^\ r=1 = 0, we may exclude the 
possibility of a saddle-node bifurcation. Furthermore, since 

Fxx(£, V) = (0, -Zim - Z1V2 + Vitif, (142) 

we may calculate the quadratic coefficient g 2 from Table [j] as 

g 2 = ±u*-F xx (u,u) = u*-(0,0, 1) T = 0, (143) 

which excludes the possibility of a transcritical bifurcation. In order to determine the third order 
coefficient 53 , we must first find the coefficient vector I120 . From Table |j] we get 

J-h 20 = -(0,0,l) T , (144) 
u*-h 20 = 0, (145) 



implying that I120 = (0, 0, |) . Observing that the Taylor expansion of (138) terminates at second 
order (F xxx = 0), we may determine the cubic coefficient 173 from Table l| as 

93 = u*-F xx (u,h 20 ) = u*-(0,-i,0) T = - b{1 S +s y (146) 
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Unfolding terms 

To find the coefficient o\ in the unfolding of the amplitude equation associated with pitchfork 
bifurcation (see Table |]) , we first observe that F r = implying that the coefficient vector hoi 
vanishes implying that the term u* F xx (u, hoi) in o\ disappears. Consequently, we get from Table |l| 



'0 0^ 

u\ = u*F xr u = u*- | 1 

,0 0, 



■u 



(147) 



For s ^ 0, we have 53 7^ , and iti ^ 0, and therefore conclude that a pitchfork bifurcation occurs 
in the Lorenz system (138) at r = 1. If s > and b > 0, the bifurcation will be supercritical. 

Fro m Ta ble [y, we now obtain the following expression for the amplitude equation for the Lorenz 
system fll3§l ) 



1 



(r - l)z - \z :: 



which admits the non-trivial stationary solutions 

z± = ±y/b(r-l). 



(148) 



(149) 



To determine the branching behavior of (|138| ) in the physical (xi, X2, o^j-space, we use the 
amplitude transformation associated with the pitchfork bifurcation from Table We find 



(150) 



We have not shown the calculation of the coefficient vectors I130 and hn, since these accidentally 
cancel when inserted into the amplitude transformation. 




A. 2 The Gray-Scott model 

As an example of a two-dimensional chemical system that exhibits a Hopf bifurcation, we con- 
sider the Gray-Scott model [||| for the transformation of a reactant A to a product B with two 
intermediate species Xi and X 2 : 

A 

Xi + 2X 2 
X 2 

Introducing dimensionless concentrations by x\ = \Jhifk\ [Xi] and xi = ^Jk^jhi [X 2 ] , a dimen- 
sionless time by r = kit, and choosing fc3 = k\ the scaled kinetic equations become 




v — x\x\ 

X\x\ — X2 



(151) 



where v = ^/fc 2 /fci[A] is used as b ifurcation parameter. Here we consider the stationary solution 
x s = (i, v). The Jacobian of (151) at x s is 



3{v) 



(152) 
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For v = 1 the Jacobian J has eigenvalues Ai = A2 = i implying that the system exhibits a Hopf 
bifurcation. The first order terms are given by the right and left complex eigenvectors u and u* 
associated with Ai. We get 



u=(-2,l + i) T and u* = + i,2i) 1 



(153) 



The second and third order terms can be obtained from Table g by calculating F xx and F x 
at the bifurcation point where v = 1 and x s = (1, 1) T . By definition 



F«(€,»j) = 



d 2 F 



dxidxj 
1,3=1 J 



tWj and Fxxx(£, tj, C) = 



d 3 F 



From (151) and (154), we get 



8 2 F 



dxf 

Second order terms 

From (154) and (155) we get 



a. 



d 2 F 



dx\dx2 



d 2 F 



8x28x1 



i,j,k=l 



8 2 F 



dxidxjdxk 



dx 2 



^■Cfc. (154) 



(155) 



* 200 = |F xx (u,u) = (8 + 2i,-8-2i) T , 
*iio = F xx (u,u) = (4,-4) T . 

The coefficients h 20 o, h 02 o, and hn can then be found by solving the linear equations 

(J-2iI)-h 200 = -*200 = -(4 + 2i,-4-2i) T , 
J-huo = -*uo = -(4,-4) T . 

giving 

2 



1200 — "020 



hno = (4,0)- 



-(5i,2^4i) J 



(156) 
(157) 

(158) 
(159) 



(160) 
(161) 



Third order terms 



From ( |15jD , (J15J) , (g_60[) , and (0_6l|) , we get 

F xx (u, h 200 ) = --(3 - 11 i)v, F xx (u, h uo ) = -(2 + 6i)v, F xx (u, h 200 ) = --(1 - 7i)v, 



F xxx (u,u,u) = 24 iv, F xxx (u,u,u) = (16 + 8i)v, 



(162) 



where v = (1,-1). The third order coefficients can then be found by solving the following linear 
equations. Solving 



(J - 3iI)-h 30 o = -$300 = -|(3 - 14i, -3 + 14i) T 



(163) 



gives 



h 30 o = h 03 o = --(-45 + 5i,42 + 9i) 
6 



(164) 
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Solving 



gives 



(J-iI)-h 210 = -* 210 = --(l-4i,-l + 4i) T , 
u*-h 2 i = 0, 



h 2 io = hi 20 = — (3 + 5 i, —4 — i). 



Using Table H, we may find the resonant coefficient 33 from (163) as 



g 3 = u* * 2 



(165) 



(166) 



(167) 



Here we notice that the resonant coefficient satisfies Re 33 < for all v showing that the Hopf 
bifurcation is always supercritical. 



Unfolding terms 

We first note that 



dF 

dv 



implying that hooi can be found from the equation 

J-hooi = -(1,0) T , 

giving 

hooi = (-1,1) T 

Observing that F XI , = 0, we may calculate the coefficient o\ as 

(Ti = u*-F xx (u, h oi) =-l + i. 
Finally, we may find the coefficient vector hioi by solving 

(J-iI>h 101 = (4,-4) T , 
u*-hioi = 0, 

which gives 

hioi = h ii = (1 - i,i) T . 

Amplitude equation 



(168) 

(169) 
(170) 
(171) 

(172) 

(173) 



In conclusion, the amplitude equation ( |100| ) associated with the supercritical Hopf bifurcation in 
the Gray-Scott model becomes 



z = iz — (1 — \)vz — (1 + - i)-z|z| 2 . 



(174) 



Any solution to the amplitude equation can then be transformed back to the concentration space 
by using 



-1 
1 

—45 + 5i 
42 + 9i 



1 + i 



5i 
2 — 4i 



z + c.c. 



\z\ 2 - 



/3 + 5 i\ . |2 , 






^_ 4 _.Jz|z| 2 + c.c. 


+ 











vz + c.c. 



(175) 



where c.c. denotes the complex conjugate of all terms within a particular square bracket. 
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A. 3 The Henon map 

As an example of an iterated map which exhibits a period doubling bifurcation, we consider the 
Henon map |34| 

Xl \ _ (a-x\ + bx 2 \ (17g) 



X 2 J \ Xi 

and derive an ex plicit amplitude equation for it using the parameter a as bifurcation parameter. 
We observe that (|l7fi| ) admits the non-trivial fixed point solution 

1 



xi = x 2 = ^{b-l+ v^a-^-l) 2 ), (177) 



with corresponding Jacobian 



J( fl )=( 1 - fe -v/4a + (6-l) 2 ) ^ (17g) 
Clearly, J (a) has an eigenvalue A = — 1 if 

The right and left eigenvectors associated with the critical eigenvalue are 

u=(-l,l) T , (180) 

(-!,&)• (181) 



1 + b 



Second and third order terms 

Observing that F xx (^,rj) = (— 2£i??i, 0) T , we find 

iF xx (u,u) = (-1,0) T . (182) 
The second order coefficient vector I120 can then be found from Table ^ by solving the linear equation 

(J-I)-h 20 = (1,0) T , (183) 

namely 

1 



hM =2(6 = i)< 1 ' 1 > ■ (184) 

Noting that F xxx = 0, we use Table || and calculate the third order resonant term g 3 as 

.93 = u*-F xx (u,h 20 ) = u*.(^-,0) T = (185) 

Unfolding terms 

To find the unfolding of the amplitude equation associated with the period doubling in ( |176] ), we 
first determine the constant coefficient vector hoi . From Table [| we get 

(J -I) -hoi = -F a = -(1,0) T (186) 
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which readily is solved for hoi as 

box = 2(^(1, If- (187) 

Here we note that all elements of the matrix F xa are zero. Using Table ||, we conclude that the 
unfolding coefficient a\ becomes 

ox = u* F xx (u, h i) = -^-j. (188) 

If we furthermore define the distance fi from the bifurcation point as /i = a — , we arrive at 

the following amplitude equation for the period doubling bifurcation in the Henon map (176) 

z ^^^—^- l ) z + Y~^ z3 - (189) 

For the doubly iterate of the amplitude equation, we find the following two non-trivial fixed point 
solutions 

z± = ±yfii. (190) 

For the Henon map, this corresponds to the emergence of a two stable period-two orbits whose 
location in (x±, £2)-phase space and parametric dependence on fx is determined via the amplitude 
transformation from Table W as 



2 

± 



(191) 



Just as in the Lorenz model (138), we have not shown the calculation of the coefficient vectors 
h3o and hn, since these accidentally cancel when inserted into the amplitude transformation. 
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Captions of tables and figures 

Figure |]: Geometrical illustration of the parameterization of the center manifold by (||) . A point 
z G E c is mapped by the transformation h to the point h(z) so that z + h(z) G W c . 

Table |l|: Expressions for calculating the coefficients of the amplitude transformation and the 
amplitude equation for each of six bifurcations: saddle-node, transcritical, and pitchfork bifurcations 
for flows and iterated maps. At the bifurcation, the Jacobian J has right and left eigenvectors u and 
u* respectively corresponding to a critical eigenvalue A = for flows and A = 1 for iterated maps. 
The amplitude transformation x = z + h(z, fj,), z = uz, transforms a solution z(t) of the amplitude 
equation to the motion x(t) on the unfolded center manifold for the dynamical system. The vector 
coefficients h pq written as h p ) are determined as solutions to the linear equations indicated, in 
terms of the derivatives of the vector field F. The coefficients of the amplitude equations can then 
be found through the explicit expressions indicated, in terms of F and h pq . 

Table ||: Formulae for calculating the coefficients of the amplitude transformation and the ampli- 
tude equation for the Hopf bifurcation and the Neimark-Sacker bifurcation. At the bifurcation, the 
Jacobian J has two complex conjugate eigenvectors u and u and left eigenvectors u* and u* corre- 
sponding to critical eigenvalues A = ± i u>o for flows (Hopf bifurcation) and A = e^ 1 16 ° for iterated 
maps (Neimark-Sacker bifurcation). The amplitude transformation x = z + h(z,yu), z = uz + uz, 
transforms a solution z(t) of the amplitude equation to the motion x(t) on the unfolded center 
manifold for the dynamical system. The vector coefficients hyfc are determined as solutions to the 
linear equations indicated, in terms of the derivatives of the vector field F. The coefficients of the 
amplitude equations can then be found through the explicit expressions indicated, in terms of F 
and hijk. 

Table [||: Expressions for calculating the coefficients of the amplitude transformation and the 
amplitude equation for the period doubling bifurcation. At the bifurcation, the Jacobian J has 
right and left eigenvectors u and u* respectively corresponding to a critical eigenvalue A = — 1 . The 
amplitude transformation z + h(z, /i), z = uz, transforms a solution z(t) to the amplitude equation 
to the motion x(f ) on the unfolded center manifold for the dynamical system. The vector coefficients 
hpq are determined as solutions to the linear equations indicated, in terms of the derivatives of the 
vector field F. The coefficients of the amplitude equation can then be found through the explicit 
expressions indicated, in terms of F and h pq . 

Table ^: Table showing the limit set solutions (stationary points, fixed points, or period orbits) 
arising from: a) saddle-node, transcritical, and pitchfork bifurcations in flows and iterated maps, b) 
Hopf bifurcations in flows, and c) period doubling and Neimark-Sacker bifurcation in iterated maps. 
The second column shows the well-known solutions to the amplitude equation, whereas column three 
gives the physical representation of these in terms of the amplitude transformation. Since the period 
doubling in an iterated map corresponds to a pitchfork bifurcation in the doubly iterated map, the 
results for these two bifurcations are identical as indicated in the table. The solution shown for the 
Hopf and Neimark-Sacker bifurcation corresponds to a periodic solution for the Hopf bifurcation. 
The limit cycle corresponds to the invariant ellipse emerging at the Neimark-Sacker bifurcation. 
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Bifu rcation 


Saddle-node 


Transcritical 


Pitchfork 


Tra nsformation 

x = z + h(z, fi) 


x = u2 + h 20 2 2 + h iM 


x = uz + h 20 2 2 + hoi^i + huz/l 


x = uz + h 2 o2 2 + h 30 2 3 + hoi/^ + huz// 


Linear equations 
for h P9 


J-h 20 = -• §Q-Fxx(u,u) 

z : 

u* h20 = 

J- h i = -Q F fl 
u*-h i = 


J-h 20 = -|QF xx (u,u) 
u*-h 2 o = 

J-hoi = -F M 
' u* h i =0 

J-h n = -Q (F XM u + F xx (u, h i)) 
' u* hn = 


J.h 20 = -iQ-Fxx(u,u) 
u* h20 = 

3 J-h 30 = -Q-F xx (u, h 20 ) + |Q-F xxx (u, u, u) 
Z : u*h 30 = 

J-hoi = -F M 
u*-h i = 

J-hn = -Q-(F XM -u + F xx (u,h i)) 
' u* hu = 


Resonant coefficients 


92 = |u*-F xx (u, u) 
(T = u*-F M 


92 = |u*-F xx (u, u) 

(Ti = u* -F XM -u + u* F xx (u, h i) 


93 = u* F xx (u, h 2 o) + jjU* F xxx (u, u, u) 
(Ti = u* F XM u + u* F xx (u, h i) 


Amplitude equation 


z = a fi + g 2 z 2 

Z CTQ/i + 2 + 


z = <j\nz + g 2 £ 2 
2h(1 + cri/i)^ + g 2 2 2 


Z = CTlflZ + g 3 z 3 

2 H-» (1 + (Tl/i)2 + g 3 2 3 


J = J for flows, 
J = J — I for iterated maps, 
Q x = x — (u* -x)u. 
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Hopf and Neimark-Sacker Bifurcations 


Transformation 

x — z -f- h(z, ^i) 


X = 


uz + uz + h2002 2 + huo \z 2 + h ( )202 2 + h3ooz 3 + h2io \z\ 2 z + 
hi2ok| 2 2 + h 302 3 + h oiM+ (hioiz + h 011 z)/i 




z 2 


(J - 72l)' 1*200 
ho20 


= -|Fxx(u,u) 

= h200 




\z\» 


(J - 7ol)-hno 


= -F xx (u,u) 




z 3 


(J - 73l)-h 3 oo 

ho30 


= -F xx (u, I1200) - |F xxx (u, u, u), 

= h300 


Linear equations 

for h ijk 


z\z\ 2 


(J - 7iI)-h 2 io 
U* -11210 
hl20 


= -Q-(F xx (u, huo) + F xx (u, h 20 o) + |F xxx (u, u, u)) 
= 
= h2io 






(J - 7oI)-h oi 


= -F M 




liz 


(J - 7iI)-hioi 
u* -hioi 

hoii 


= -Q-(F XM -u + F xx (u,h oi)) 
= 
= hioi 


Resonant coefficients 


93 -- 

(71 = 


= u* • F xx (u, huo ) + u* • F xx (u, I1200 ) + \ u* ■ F xxx (u, u, u) 
= u* F x(J u + u* F xx (u, hooi) 


Amplitude equation 


2 = 
Z 1— > 


(iu! + (J\p.)z + g 3 z\z\ 2 
(c ie o +a lf i)z + gzz\z\ 2 


7fc = fciuo for flows, 
7fe = e kl6 ° for iterated maps, 
Qx = x — (u* x)u. 
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Period Doubling Bifurcation 


Transformation 

x = z + h(z, fi) 


x=u2 + I120Z 2 + h3oz 3 + hoiA* + hn^ 


Linear equations 
for h P9 


2 2. (J-I)-h20=-|Fxx(u,u) 

(J + I)-h 30 = -Q-(F xx (u, h 20 ) + |F xxx (u, u, u)) 

z ■ 

u -h3o = 
(J — I) -hoi = — F M 

(J + I)-hn = -Q-(F XM -u + F xx (u,h i)) 
u*hn = 


Resonant coefficients 


93 = u* F xx (u, h 2 o) + |u* F xxx (u, u, u) 
cri = u* -F XM -u + u* F xx (u, hoi) 


Amplitude equation 


z i-> (ai/j, - l)z + 93Z 3 


Q x = x — (u* x)u. 
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Limit set solutions 


Bifurcation 


Amplitude space 


Physical space 


Solution diagram 


Saddle-node 




x± = ±ud + h 2 oa 2 + h i/x 






1 I 


z± = ±a, a = . / 

V 92 






1 S\ 








"t 




Transcritical 


z + = a, a = 

52 


x + = ua + h 20 a 2 + h i^ + hna/i 
x = h i/z 




Z 


/\ 


(..•-" 

-"i 


1 
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Period doubling 




x± = ±ua + h 2 oa 2 ± h 30 a 3 + h i/x ± hya/i 
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